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Abstract 

Lifespan distributions of populations of quite diverse species such as humans and 
yeast seem to surprisingly well follow the same empirical Gompertz-Makeham 
law, which basically predicts an exponential increase of mortality rate with age. 
This empirical law can for example be grounded in reliability theory when in- 
dividuals age through the random failure of a number of redundant essential 
functional units. However, ageing and subsequent death can also be caused by 
the accumulation of "ageing factors" , for example noxious metabolic end prod- 
ucts or genetic anomalities, such as self-replicating extra-chromosomal DNA in 
yeast. 

We first show how Gompertz-Makeham behaviour arises when ageing factor 
accumulation follows a deterministic self-reinforcing process. We go then on to 
demonstrate that such a deterministic process is a good approximation of the 
underlying stochastic accumulation of ageing factors where the stochastic model 
can also account for old-age levelling off of mortality rate. 
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1. Introduction 

Very simple organisms do not age, that is they do not show a deterioration of 
significant life functions with age. Instead they seem to simply break, when one 
of their essential functions happens to fail. If such failures occur independently 
with a constant rate A - that is failures follow a Poisson process - then their cor- 
responding survival function S (t) decays exponentially with time, and lifespans 
are distributed exponentially - this is just a simple conclusion from mathemat- 
ical reliability theory and other formally similar processes such as radioactive 
decay. 
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More complex organisms however do age, ie show a deterioration of hfe func- 
tions, for example regarding regeneration after injury, reproduction, metabolic 
activity and so on, and their mortality rates X{t) tend to increase in time. For 
many different species, it has been empirically established that mortality rates 
X{t) increase exponentially with time, following the famous empirical Gompertz- 
Makeham law of mortality 

A(t) = Ao + ae''*,Ao,a,^ > 0, (1) 



at least approximately (jGavrilov and Gavrilovalll99lh where a,/? are constants 



and Ao the age-independent death rate. Notable deviations are usually for very 
young ages t when "childhood diseases" lead to higher mortality and also for old 
ages where empirically the mortality rate increases slower than exponentially or 
even levels off to a constant value. 

In this paper we will assume that the driver of ageing, and ultimately death, 
is the accumulation of Age ing Factors (AFs), wh ich for example plays a major 



role in the ageing of yeast ([Sinclair et al.l . Il998al ). We will demonstrate generi 



cally that when such AFs are produced in an auto-catalytic or self-reinforcing 
process, then Gompertz-Makeham behaviour of the mortality rate follows. Such 
self-reinforcing processes are quite generic, for example in auto-catalytic chemi- 
cal or biological reactions or also in self-replicating genetic anomalities, and show 
typically exponential behaviour. The term Ageing Factor (AF) will therefore 
generically stand for any substance or anomaly that is thought to cause ageing 
in a specific organism. AFs could come in discrete (even "macroscopic" ) chunks 
or be a continuous quantity. They could be measured in absolute numbers of 
Ageing Factor Units (AFU s) or as a concentration. 

We proceed as follows: We will firstly take the naive point of view that 
creation and replication of AFs is a deterministic process and mortality propor- 
tional to the abundance of AF. This straight-forwardly leads to an exponentially 
increasing mortality rate (with no old-age levelling). We then go an and assume 
that the creation and replication of the AFs are stochastic when the AF comes in 
discrete units and in small numbers. If mortality rate is assumed to be propor- 
tional to the probability that an organism has an AF abundance above a critical 
level, then it follows that the mortality rate increases (more than) exponential 
for middle ages while it levels off to a constant value at old-ages. 



2. Background 

2.1. Mortality Models 

The above equation eq. ([T]) is primarily descriptive and empirically derived 
from data, however does not in itself indicate a reason why mortality should 
increase exponentially with age (at least for middle ages) or any underlying 
mechanism(s) of ageing. However any such mechanism would need to be quite 
general since lifespan of so many different species often surprisingly well fit the 
Gompertz-Makeham law. Therefore various attempts have been undertaken to 
ground Gompertz-Makeham behaviour in underlying processes that could give 
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an explanation as to why the mortality rate should increase exponentially with 
time. 

One prominent such attempt follows IGavrilov and Gavrilova (|200lh . These 



authors ground the Gompertz-Makeham (and Weibull) laws in reliability theory. 
Their basic assumption is that complex organisms have a block of redundant 
functional elements for each essential function. Each such element does not age 
and hence fails with a constant rate. The organism as a whole fails (or dies) if 
all redundant elements for a single function have been exhausted. Depending on 
the number of redundant elements and their failure state at birth, a Gompertz- 
Makeham or Weibull law of mortality results. This model of ageing also nicely 
explains the levelling off of the mortality rate at old ages. 

A second attempt to grou nd the Gompertz-Makeham law in a biological 



mechanism is lShklovskiil (|2005f) starting from the assumption that the ability of 
on organism to neutralise defective cells (or noxious substances for that matter) 
is based on random encounters of these harmful items with some neutralising 
antagonist. If the expected number of such encounters sinks linearly with time, 
again the Gompertz-Makeham law follows. However no attempt is made to 
explain old-age longevity or to make plausible why the frequency of such random 
encounters should decrease linearly in time. 

While the two models above are individual-based, ie when an individ- 
ual dies is determ i ned by its l i fe his t ory, the r e is f inally the Penna model 



( Laszkiewicz et all 2005 . Penna . 19951 Stauffei . 2007t) . This, too, can explain 



an exponentially increasing mortality rate, albeit on a population level. The 
Penna model assumes that genetic defects come to bear at different ages, and 
individuals die if they have accumulated a certain number of genetic defects - 
in that sense it is similar to ageing in yeast below. However how many defects 
have been acquired and when they come to bear is predetermined genetically for 
each individual, and only a population equilibrium of individuals with different 
defective genes leads then to a Gompertzian exponentially increasing mortality 
rate. 

2.2. Ageing in Yeast 

Our approach is motivated by systems that age through accumulating an 
AF. Yeast (Saccharomyces cerevesiae) has in this respe ct been studied exten - 



sively and served as model for ageing in higher organisms (jSinclair et al.l . ll998al ). 
Albeit a single-cell organism, yeast divides asymmetrically, so that mother and 
daughter cells can be identified. Mothers eventually cease producing daughter 
cells, stall and subsequently die. Extensive experimentation has shown that 
so-c alled Extra chromosom a l rPN A circles (ERCs) are a probable cause of age- 



ing ( Sinclair and Guarentd . [T997h . Finally, lifespan of yeast is best measured 



in number of generations, that is how many times a mother cell has already 
divided. While the absolute time between two such divisions depends on the 
environmental conditions, lifespan distribution measured in generational age is 
constant for different environments. 

An ERG is a snippet of a repetitive section of chromosomal 
rDeoxyribonucleic Acid (DNA) that is excised from the chromosome sponta- 
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neously with a low rate, and subsequently forms an extra-chromosomal DNA 
ring. Excision can happen several times independently with approximately con- 
stant low rate since the DNA section in question contains about 100 repeats 
of the snippet. As the snippet has its own begin-of-replication sequence, once 
excised it replicates with a certain probability in each cell cycle - and hence 
replicates in a self-reinforcing process. 

At cell division, the mother cell retains almost all of the ERCs and does 
not share them - eg proportional to cell volumes - with the daughter cells, so 
that ERCs accumulate in mother cells whereas daughters usually start life free of 
ERCs. Cell death then is assumed to occur because replication of a high number 
of ERCs (in the order of 1000) exhausts cell resources needed for replication of 
the core chromosomal DNA. 

For very old mother cells ERC retention capability seems to saturate, so 
that ERCs are shared with daughters, effectively decreasing the mother's ERC 
replication rate and producing prematurely aged daughters with a reduced 
life span. Modelling the effects of retention saturation is beyond the scope 



of this paper. From figures discussed in the liter ature (jGillespie et al.l . 12004 . 



Sinclair and Guarentd . 119971 [Sinclair et al.l Il998bl ). it seems however that the 



reduction in ERC number is minor compared to their absolute number and/or 
rate of replication, so that mothers' lifespans are extended around one or two 
generations at most. It does however significantly reduce the expected lifespan 
of daughters of old mothers that have received a portion of the mother's ERCs. 
The model discussed below can straight-forwardly be modified to take into ac- 
count birth with a number of ERC but for sake of simplicity we will only deal 
with populations that start life from a clean ERC- fr ee state as has been the case 
in many experiments (jSinclair and Guarentd . 119971 ) . 

Experiments suggest as well that for yeast, cells live on quite happily when 
they have low or medium numbers of ERCs but die when they have reached a 
critical level of ERCs so that there does not seem to be a linear increase of mor- 
tality rate with number of ERCs as we will naively assume in the deterministic 
model. Instead the underlying assumption seems to be that mortality s hows an 



increase at or around a critical number of ERCs (jSinclair et al. 

To conclude, agein g (in yeast and 

Sinclair. Mills, and Guarentd Il998al) can be caused by 



1998bD . 



other 
an AF 



organisms, 
that comes 

in discrete chunks and where individuals initially start up with no or a very 
small number of such AF. The number of such AF increases due to two 
processes: initial creation of such AFs (in yeast by excision from chromosomal 
DNA) and subsequently self-reinforcing replication (for yeast once in the cell 
cycle). Finally mortality is correlated with the abundance of such AFs. 

In the next section we will present two formalisations of these processes. The 
first is naive, straight-forward and deterministic. The second however is more 
complicated and based on stochastic processes. This stochastic model serves 
also to justify the naive approach. These approaches are certainly motivated by 
the experimental data for yeast, but we will formulate them quite generally as 
auto-catalytic processes are an ubiquitous phenomenon in nature, thus might 
well play a role in ageing processes in many organisms. 
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3. Models 

Let c{t) denote the abundance of an AF at time t. The abundance could 
either stand for the (discrete) absolute number of Ageing Factor Units (AFU 
s) if the AF comes in discrete chunks and absolute numbers are low, or for 
a (continuous) concentration if absolute numbers are high or the factor itself 
is a continuous variable. We shall assume that c(t) changes in time due to 
two processes: 1. Creation with rate p(t), ie new AFU s come into existence 
independently of existing ones. This could be from internal processes such as a 
product of a chemical or metabolic reaction or the excision of snippets of DNA 
as in yeast, but also external such as a DNA defects due to UV light and so 
on. 2. Replication with rate r{t), ie once some AFU s have been acquired there 
is an auto-catalytic process that produces more AFU s, linearly depending on 
the abundance c{t) of AF already present, ie Ac{t) ^ c{t). Typical examples of 
such processes are auto-catalytic chemical reactions where the product at the 
same time is also an educt of the reaction, or the replication of a DNA snippet 
independently of the core genome as happens with the ERCs in yeast. These 
processes could be happening in continuous or discrete time, on an absolute time 
scale or tied to an internal time scale. For yeast for example, ERC replication 
is synchronised with the cell cycle. 



3.1. Deterministic Production of AFs 

As it is simpler and instructive, we first consider the case where the AF comes 
in large quantities in each individual and creation and replication processes have 
high rates. This means that while individual creation and replication may well 
be stochastic, we can assume - due to the high rates and high numbers of AF 
involved - that creation and replication proceed deterministically, and variations 
of AF around its instantaneous mean are negligible. 

We will in addition assume continuous time. Mathematical argument is often 
easier in continuous time than in discrete time, and the inaccuracies incurred 
by transforming equations between discrete time and continuous time are small 
when the involved quantities are smooth. 

Formalising creation and replication processes in continuous time, c{t) fol- 
lows an inhomogeneous linear Ordindary Differential Equation (ODE): 

^^r{t)c{t)+p{t) (2) 
with solution for initial value c(0) = cq (see [Appendix A ): 

c{t)^ f p{T)e^>^^'^'^''dT + coeIor^-}'ir 

c{t),p{t),r{t) are non-negative functions, and it is obvious that c(t) increases 
essentially exponentially with t when r{t) > 0. For discrete time, c(t) would be 
the solution of a similar difference equation which again increases exponentially 
(geometrically) in time. 
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For simplicity let us assume r{t) — r and p{t) = p are constant, so that both 
creation and replication processes of AF do not change with time, and eq. ([3]) 
simplifies to 

c(t) = ^(e'-*-l) + coe'-* (4) 

However time-varying r(t) and p{t) could also easily be accommodated and 
would not essentially change the arguments belows. 

If we assume the death is only due to the AF, then a natural assumption is 
that the mortality rate is proportional to c(<), ie: 

X{t) ^ zc{t),z>0 (5) 

where z is the proportionality factor. Then quite straight-forwardly, inserting 
eq. (U]) into eq. ([5]) provides a Gompertzian time development of t with some 
additional terms to adapt to initial conditions. For higher t, X(t) increases 
exponentially with t, and does not level off to a constant value either for very 
high t. 

3.2. Example: Lifespan of Yeast 

Experimental data is often in the form of lifespan distributions rather than 
mortality rate per se. The lifespan distribution is given by the survival function 
S{t) which is the probability that an individual survives at least until t. It is 
related to the mortality rate \{t) via the following ODE with initial condition 
^(O) = 1: 

dS{t) 



dt 



-\{t)S{t). (6) 



This is again of the form as in Appendix A| and inserting eq. ([4]) and eq. ([5]) 



and solving for S{t) (see eq. (|B.2|) ') we get: 

S{t) ^ e-SoHr}dr ^ (7) 

This form of S{t) can now straight-forwardly be fitted to experimental data, see 

fig.m 

As p and cq do not enter eq. (O independently of z, their absolute value 
cannot be estimated from data unless further biological assumptions are taken 
into account, r can be seen as the average or effective rate of increase of the 
AF. The fit is surprisingly good even though it does not take into account the 
levelling of the mortality rate X{t) observed for many species for high t so it 
underestimates the frequency of very high lifespans. 

3.3. Stochastic Production of AFs 

While eq. ([7]) often yields a good fit as we have seen, it starts from the 
assumption that the AF comes in high numbers with high rates p and r, so 
that the process of AF is largely deterministic. The deterministic model does 
also not account for the levelling off of mortality rate with time. For example 
for yeast it is however known that the AF comes only in discrete units and 
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Experiment 
Deterministic IVIodei 



Figure 1: Fitting the surviva l curve of yeast: A least squares fit for p' := zp and r of eq. (0 
against [Sinclair etaP s data Il998bl (Fig. 1). CO was set to 0, as experiments started with 
ERC clean cells. Also no age-independent mortality rate Aq was used, as accidental death 
is often negligible under lab-conditions. Best least squares fit for p' = 0.00011, r = 0.22, 
SSR = 0.0034, RMS = 0.014. Note the (expected) underestimation of old-age longevity. 



individuals acquire them one by one, and numbers stay smallish throughout 
life-time. All this does not warrant a deterministic approach. Hence for such 
discrete AFs in small numbers, creation and replication should be treated as 
stochastic processes. For simplicity of argument here, we will assume discrete 
time in the following. 

The stochastic creation process of discrete AFU s is easily modelled as a 
(constant or varying rate) Poisson or Bernoulli process, so for constant rate p 
the expected number of AFU increases as pt with variance pt. In the following 
we will assume that p is constant. This is simpler and there is also no biological 
indication to contrary. 

The replication process is more complicated and be st be described in the 



framework of branching processes (|Haccou et al.l . 120051) . a type of stochastic 



process suited to model growth processes. We will first look at a pure branching 
process that corresponds to the stochastic AFU replication in isolation, and 
subsequently combine replication and creation processes. 

In the following capital letters such as X denote random variables. Pr(-) 
denotes the probability (density) of an event, E(-) and Var(-) the expectation 
value and variance of a random variable respectively. For further background 
- such as relations of conditional expectations and variances or approximations 
of the Q- and erf-function - we refer the read er to good textbook s with an 
introduction t o pro bability theory, for example ICover and ThomasI (jl991r ) or 
Haccou et all (|2005l ). 
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3.4- Replication Process Only 

In this section we concentrate on the rephcation process per se, and assume 
that initiahy one AFU is already present: c(0) = 1. Then existing AFU s 
rephcate with mean rate r in each time step. This means each single AFU 
at t gives rise to one (itself) or two AFU s (itself and its replication) at t + 1 
with probabilities 1 — r and r respectively. Formally individual replication is a 
random variable S with probabilities Pr(S = 1) = 1 — ?- and Pr(S ~ 2) ~ r, and 
all other outcomes have probability zero, and we define short symbols for mean 
and variance of the number of successors of a single AFU : m := E(S) — 1 + r 
and :— VarCE.) — r(l — r). 

As the replication process itself is random, also the abundance of the pop- 
ulation c{t) at time t is a random variable and - to mark that the abundance 
is now a random variable in discrete time - we replace symbol c(t) with Ct and 
call it the population size. 

In fact, the sequence of (Ct) forms a branching process (see [Appendix C ) 



and Ct is the Ct_ i-fold sum of the individual replication process :z (with i.i.d. 
outcomes; we refrain from distinguishing the individual S in order not to over- 
burden the notation): Ct = X]i=i^ ^- The expectation value /i(i) :— E(^) and 
variance cr^(f) := Var{Ct) of the population sizes Ct then follow from expecta- 
tion value and variance of individual offspring S as (see [Appendix C ): 



2 m^{vn} — 1) 



M(t) = m*, a\t) = s^^ (8) 

m[m — 1) 

with the initial condition that we start from a single AFU at < = 0, ie Co = 1. 
In the case of deterministic growth, the mortality rate \{t) was assumed 



proportional to cit). We could in principle also here set \{t) ^ E(Ct) = m* 
with the same results as above, loosing again the levelling off of \{t) for old 
ages. Experiments indicate that for yeast, cells with a relative low number of 
ERCs do not have an altered instantaneous mortality rate from cells with no 
ERCs. The mortality rate seems to increase only when the number of ERCs 
has reached a critical level. Hence in this stochastic process setting, it is more 
natural to assume that the mortality rate X{t) is proportional to the probability 
that an individual has more AFU s Ct than a fixed critical level c-f , ie 

A(i) = zPr(Ct > ct),z > 0. (9) 

This agrees also wit h the type of mortality criteria used in computer simulations 

ancies are 



i nis agrees also wit n tne type or mortanty criteria usea m computer sim 
of ageing processes ( Gillespie et al. I . [20041 ). More complicated depender 



of course possible, but would typically show a similar dependency on an upper 
quantile of Ct- 

In order to calculate \{t) explicitly, we would need to know the distribution of 
Ct for all times t. To our knowledge this problem has not been solved explicitly. 
We can however make a series of approximations for small times t, large times t 
and middle times t to demonstrate that this setting for the mortality rate X(t) 
yields Gompertz-Makeham behaviour for middle t and a constant mortality for 
very old ages. 
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Small t. We assume that at time t = there is one AFU , Cq — 1. As each 
AFU can have at most two successors (itself and its offspring), we know that 
for t < log{c-\)/log{2), Pr(Ct > cq) = 0. So depending on there wiU be more 
or less long an initial period for which the mortality rate is zero. 

Approximation with normal distribution. Generally the precise form of Pr (Ct) 
is not known. Knowing only its /i(i) and a(t) we assume that Ct follows ap- 
proximately a normal distribution with mean /i(t) and variance cr'^it). This is 
plausible as Ct ultimately in our setting is a sum of binomial distributions (each 
individual either replicates or not and survives itself into the next generation, 
so Ct+i = Ct + B{Ct,r)) which themselves can be approximated by a normal 
distribution - however the Ct are not independent. Also quantities related to 
Ct such as Pr(Ct|(7t i) hav e been shown to approximate a normal distribution 
( Quine and Szczotka . Il993). Hence for mildly large t we approximate: 



p.(c,.c,,»<,(a_M).^4e.(a^) ao, 

where Q{x) — 1/2— 1/2 erf(a;/-y/(2)) is the Q-function that gives the mass of 
a standard normal distribution above x (ie the probability of getting a value 
larger than x). 

Large t ^ 0. For large very large t ^ 0, 

^.^ c^-MO ^ ^.^ cy^' ^ _ 1 

t->oo a{t) t->oo / in*(m*-l) S 

Y m(m — 1) 

as both ^{t),a{t) grow exponentially with the same rate. Hence 

lim X{t) = Qi-l/s) > 0, (12) 

ie the mortality rate X{t) levels off with age to a constant positive value, hence 
X{t) shows the experimentally observed old-age deviation from the Gompertz- 
Makeham law. 

Middle t. We have established that X{t) — for small t and w const > for 
large t. As each AFU survives into the next time step, Ct+i > Ct, ie the 
number of AFU s is non-decreasing. This means that X{t) — zPr{Ct > c-f) is 
non-decreasing as well. Hence X{t) has roughly the form of a (discrete) sigmoid. 

If approximation of Ct with a normal distribution is valid also for middle 
t such that ^ fi{t) < C|, where P{Ct > C|) > but the expected number 
fj,{t) — TO* still less than the critical value C|, we go on to show that increase of 
X{t) must be steeper than to*, in other words we have an exponential growth of 
the mortality rate for intermediate ages. 

If X :— — /i(i)/(T is not too close to zero, ie /i(t) still much smaller than 
Ct, then Q(x) w — e~^^ for x > and hence 

m . g (^-1^) . ^ e-^. (13) 

^\ a{t) J V2^(ct - /.(t)) ^ ' 
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For C| > the exponential is increasing towards 1. The denominator of 

the pre-factor is decreasing, and its numerator increasing exponentially approx- 
imately as TO* (since w to* for not too small t). Hence there are 
intermediate t for which X(t) approximately increases more than exponentially: 

A(t) > zm* (14) 

3.5. Creation and Replication 

Unless born out of very old mother cells, new-born yeast cell usually do not 
have any ERC, ie Ct = so creation and replication process need to interact to 
cause ageing. Before replication of AFU s can start, first some must be produced, 
ie we have a combination of two stochastic processes. This combination is more 
cumbersome to deal with, but we demonstrate below that (asymptotic) growth 
rates of E(Ct) and Var(Ct) are unaltered, so that the same approximations and 
asymptotic considerations are valid for this more complicated process as were 
in the case of replication only. 

Let Ct denote the AFU population size as before, and S again the individual 
AFU replication process with the same distribution. Finally let Xt be the 
creation process with constant rate p, ie Pr(Xf = 0) = 1 —p and Pi{Xt — 1) ~ p. 
Then E(X) = p and Var(X) = p{l - p). 

The random number Ct+i of AFU at time t -I- 1 is then the sum of the 
replication of the previous number of AFU s Yliili " P^^^ the number Xt of 
newly created AFU s: 

Ct 

Ct+i =Y,~+Xt (15) 

4=1 

Utilising independence of the individual S and Xt we can calculate a recursive 
equation for ]E(Ct) as in the case above, see [Appendix D[ 

E(Ct+i) =toE(CO+E(XO (16) 

with explicit solution for initial value Co = 0: 

E(a) = ^(m*-1). (17) 
to — 1 

Comparison of pB]) and fTT]) with ^ and ^ also makes it clear that the process 
Ct is the stochastic time-discrete analogue to the deterministic continuous time 
equation in section [3T] (with r — log (to)). 

Likewise a recursive formula for Var(Ct) can be obtained: 

Var(Ct+i) = TO^ Var(Ct) + E(Ct) + Var(Xt). (18) 

This has a rather tedious explicit solution, see [Appendix D[ For large t however 
it is clear from the dependency of Var(Ct) on Var(Ct_i) and because E(Ct) only 
grows like to* for large i, that the growth of Var(Ct) will asymptotically be like 
TO^*. Ie both Var(Ct) and E(Ct) grow asymptotically as in the case before. 
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Hence again under the assumption that Ct can be approximated well enough 
with a normal distribution, the same considerations as above for small, middle 
and large follow, ie A(<) = for small t, X> zm* for middle t and X{t) constant 
for large t. 



4. Discussion 

The production processes and death criteria considered in this paper arise 
from experimental observation of ageing in yeast mother cells where ERCs seem 
to be drivers of ageing and death. It is remarkable indeed that the deter- 
ministic equation eq. (O already yields such a good fit with experimental data, 
even though the underlying process model is deterministic, uses a non-biological 
death criterion and finally does not take any detailed properties of the under- 
lying AF production processes into account such as sharing of ERCs between 
old mothers and daughters. It only relies on the exponential accumulation of 
AFs. The model does not take into account either that ageing in yeast need not 
necessarily only take place through accumulation of ERC. 

In the stochastic model, we have demonstrated that the mortality rate in- 
creases with time roughly like a discrete sigmoid and X{t) > zm* — ze^^^r = 
log{m) for middle t. The good fit of an exponential mortality rate with data 
(eg in fig. [J) shows that usually \{t) does not increase significantly faster than 
exponentially or that such faster increase has no significant effect on the form 
of the survival function. If we do not care about the overestimation of mortality 
old ages, we can extend the exponential approximation of A(t) also to greater 
t. Finally we know that Pr(Ct > Cf ) =0 strictly for young ages. Extending the 
approximation of Ct with a normal distribution also to these young ages would 
lead to a Pr(Cf > C|) > 0. However as then also c-f ^ m{t) this is often negli- 
gibly small (and would for example be masked by a constant age-independent 
death rate Aq). Hence the deterministic approach \{t) ^ e''* can be justified for 
all ages t as an approxim ation of the stochastic approach \{t) ^ Pr(Cf > Cj). 



Gillespie et al.l ()2004() set out to understand ageing in yeast with the help of 
numeric computer simulations. Essentially they simulated a stochastic process 
such as in section 13.51 Their simulation also modelled the old-age ERC sharing 
between mothers and daughters. Simulation parameters, such as replication 
rate r or excision rate p were estimated from data in the literature. With these 
estimated parameters, their simulation yielded a fairly good fit which however 
was not quantified. The best fit was achieved by assuming that p{t) ~ t^, ie 
a non-constant creation rate which seems to be an ad-hoc choice to improve 
the fit over a constant excision rate. In light of the go od fit of eq. ([7ll which 
was derived with a constant p, it is quite surprising that iGillespie et al.l did not 
obtain a good fit under the same assumption. 

Finally the rat e with which ERCs in yeast a re replicated has been estimated 



as about m = 1.6 (jSinclair and Guarentelll997l) - this follows as an average rate 



from an average lifespan of 15 generations once the first ERC has been acquired 
and a critical level in the order of 1000, as to « 1000^/^^ = 1.585 and hence 
r — log 1.59 sa 0.46, this is about twice the r estimated from the data in fig. [TJ 



11 



Last but not least, eq. (O and its Rood fit to the data also confirm quan- 
titatively a qualitative conjecture of Sinclair et al. ( 1998bf ), namely that the 



survival function for yeast is a "sum" of a decay process (due to excision) and 
subsequent exponential accumulation of ERCs, albeit t he combination of these 
two process is a more involved sum than [Sinclair et "aL implied, namely that of 
two stochastic processes in eq. p5|) . 

5. Conclusion 

In this paper we have demonstrated that mortality rates similar to those in 
the Gompertz-Makeham law can be derived when we assume an AF is produced 
in an auto-catalytic process. Such processes in turn lead to an exponential (or 
synonymously geometric) increase of the abundance of AFs. In the case of 
deterministic production of AFs, with mortality assumed proportional to the 
AF abundance a plain Gompertz-Makeham law follows with the usual under- 
estimation of longevity for old ages. In the stochastic case, not only does a 
Gompertz-Makeham law follow, but also the levelling off of mortality rates is 
derived if mortality is assumed proportional to the probability that an individual 
has reached a critical level of AFs. 

All in all the paper establishes a useful and interesting relation between 
mortality rates that follow a Gompertz-Makeham law and an underlying auto- 
catalytic process. For yeast ageing driven by ERCs the model and data are in 
qualitative agreement. However discrepancies between the rate of replication 
estimated here and within other experimental contexts need to be examined. 

The derived relation is general and applicable to all organisms where ageing 
is based on exp onential accumulat ion of an A F* - th is nicely contrasts and com- 
plements Gavrilo v and Gavrilovaf s approach (I2OOII) for which ageing is based 



on the exhaustion of parallel redundant functional elements. 
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Appendix A. General Linear ODE 

In the present paper the linear inhomogeneous ODE plays are role: 

dx{t) 



dt 



^ a{t)x{t) + b{t) (A.l) 



where a(t) can be thought of as a time-variant rate of growth of x{t) and b{t) 
as a time- variant external influx to x{t). The generic solution can be obtained 
- under some smoothness assumptions - for example using the Green's function 
approach. We refer the reader to a generic textbook on ODEs and just state 
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the generic solution below. The reader can check the correctness of the solution 
through differentiating with respect to t and inserting into eq. (lA.ll) . 

In the homogeneous case, i.e. b{t) = 0, the general solution of eq. (|A.1I) with 
initial value a;(0) = xq is 

xit) = xoe^o . (A.2) 

In the inhomogeneous case b{t) ^ 0, the general solution with initial value 
x(t) = Xq is: 

X{t)^ p{T)e^>'-^'^'^^'dT + Xoe^o^'^-)dr (^3) 

Jo 

For example if a{t) — a, b{t) ~ b are constant and x{0) = the solution is: 

a;(t) = -(e°*-l) (A.4) 
a 



Appendix B. Lifespans and mortality rate 

Let X{t) be the instantaneous mortality (or failure) rate, ie the rate with 
which an individual dies/ fails at t, then is it related to the distribution of lifes- 
pans or survival function S{t) (ie the probability for an individual to survive 
until time t) via: 

^ = -ms{t) (B.i) 

with initial condition S'(O) — 1. This equation is of the form of eq. (jA.l[) . hence 
solved for S{t) wc obtain: 

5(t) = e^-^'o^f^)''^. (B.2) 

If we have independent processes that cause death with rates Xi{t) and 
X2{t) and corresponding survival functions Si{t),S2{t) and combine them as 
X{t) — ciXi(t) + 02X2(1), then due to linearity of integration and because the 
exponential transfers sums into products: 

5(t) = 5JHi)^rW- (B.3) 

This is handy when dealing with different terms in a complex X{t) separately. 
For example adding constant age independent mortality rate Ao to a X{t) means 
that the corresponding S{t) acquires an additional factor e^'^°*. 



Appendix C. Branching Processes 

Intuitively a branching processes is a model for stochastic reproduction of 
individuals in discrete time and follows the population size Ct in time. Each 
individual lives for one time step (or generation), and produces a number of 
offspring according to a random variable (rv) S with mean of m :— E(S) and 
variance := Var(S) according to a fixed probability distribution. Given the 
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population size Ct, the number of individuals in the next generation Ct+i is 
a random number that is the sum of all offsprings produced randomly and 
independently by all individuals in Ct, ie Ct+i = Yl'i=Q^- The sequence of 
(Ct)t>o f orms a stochastic pr ocess, ie a sequence of interdependent random 
variables ( Haccou et al. . 20051) . 



Let E(X) denote the expectation value of a random variable X and Pr(X) 
and Var(X) its probability and variance. For given expectation m of an in- 
dividual's offspring and initial population Cq = 1, one can then calculate 
the expectation values ii{t) := E(Ct) of expected population sizes. Since 
E(Ct+i|Ct) = ES,;IE(S) = mCt it follows: 



iiit + 1) = nct+i) - nnct+i\Ct)) = 

E(mCt) =™IE(Ct_i). (C.l) 

Utilising this recursiveness, it follows with initial value Cq = 1: 

H{t) = m*. (C.2) 

A similar but more le ngthy argument yields (see for example 
Haccou. Jagers. and VatutinI [ioOSi 'l or compare to the more complex ar- 
gument below for eq. (|D.8P ) 

a\t) Var(Q) = - 1) (^.3) 

mym — 1) 

For ^ 3> 0, this can be approximated as 

a^{t) « s^m^* (C.4) 

Hence for a branching process expectation and its standard deviation a{t) 
both increase exponentially with to*. 

"Individuals" and "offspring" in a branching process as above can correspond 
to "individuals" in the biological sense but also to other discrete units that 
reproduce. In the main text, an individual is an AFU that reduplicates with 
probability r in a time step and also survives into the next time step with 
certainty. An individual AFU then has either one (itself, no replication) or two 
successors (itself and its replication), ie individual reproduction S into the next 
time step is governed by the following probability distribution: 

Pr(S 1) = (1 -r), Pr(S = 2)==r (C.5) 

and probabilities for all other outcomes zero. The individual reproduction pro- 
cess S then has mean m — 1 + r and variance s"^ — r{l — r). Note that this 
branching process cannot die out (Pr(Ci > 0,Vt) = 1), as there is at least one 
successor for each individual (ie Pr(S > 1) = 1). In fact the process will be 
non-decreasing Ct+i > Ct- 
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Appendix D. Combined Production Process 

The full stochastic production process of AFU s is the combination of a 
standard Galton- Watson branching process that governs AFU replication and 
a Bernoulli process (a discrete time Poisson process) that governs AFU creation. 

Let Ct denote the population size at i, and S the individual AFU replication 
process with mean 1E(S) — m and variance Var(S) — as before. Let further 
Xt be the Bernoulli creation process with constant rate p, ie Pr(At = 0) = 1—p 
and Pr(At = 1) = p. Then E(At) = p and Var(At) = p(l — p) at each time 
step. The process we are interested in is given by the sequence of recursively 
defined random variables Ct 

Ct 

Ct+i=^S + At. (D.l) 

1=1 

Utilising independence of the individual Ss and Xt we can calculate a recursive 
equation for E(Ct). First we note that 

Ct 

EiCt+i\Ct) = 5]E(S) + E{Xt) = Ctm + p. (D.2) 

Then making use of E{A) = E{E{A\B) it follows: 

E(Ct+i) = E(E(Ct+i|Ct) = mE(Ct) +p. (D.3) 
With initial condition Cq — E(Co) = this can be resolved in closed form: 

E{Ct) = ^-(m* - 1). (D.4) 

m — 1 

Consequentially E(Ct) increases with m* for large t. For the variance, we set 
Zt := X]i=ii hence Ct+i = Zt + Xt and utilising again the independence of Zt 
and Xt: 

Var(Ct+i) = Var(Zt) + Var(At) (D.5) 

and further using the variance partitioning formula Var(A) = E{Va.r{A\B)) + 
Var(E(A|B): 

Var(Zt) - E(Var(Zt|Ct)) + Var(E(Zt|Ct)) (D.6) 
= E{Cts^) + Var(mCt) = E{Ct) + Var(Ct), (D.7) 

and hence we get the recursive formula: 

Var(Ct+i) = Var(Ct) + E{Ct) + Var(At). (D.8) 

We have already a closed solution for E(Ct). Inserting eq. (ITTl) into (jD.8[) we 
(or rather maxima with some subsequent manual simplification) can solve this 
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for Var(Cf) with initial condition Var(Co) = (from Co = 0): 



Var(Ct) -m2*(p(l-p) 



-2t 



1 




l-m"* 




) (D-9) 



m{m — 1) 



(m — l){m? — 1) 



This is rather cumbersome, and we refrain from further simphfication as it is 
obvious that for t — >• cx) Var(Cf) ^ m^*. As also this Ct is non-decreasing, 
the same argumentation as in section 13.41 can be applied to the full process in 
section 13.51 
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